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Abstract 

We introduce SuperLFV, a numerical tool for calculating low-energy observables that 
exhibit charged lepton flavor violation (LFV) in the context of the minimal supersymmetric 
standard model (MSSM). As the Large Hadron Collider and MEG, a dedicated /i"*" — > 
e^7 experiment, are presently acquiring data, there is need for tools that provide rapid 
discrimination of models that exhibit LFV. SuperLFV accepts a spectrum file compliant 
with the SUSY Les Houches Accord (SLHA), containing the MSSM couplings and masses 
with complex phases at the supersymmetry breaking scale. In this manner, SuperLFV is 
compatible with but divorced from existing SLHA spectrum calculators that provides the 
low energy spectrum. Hence, input spectra are not confined to the LFV sources provided by 
established SLHA spectrum calculators. Input spectra may be generated by personal code 
or by hand, allowing for arbitrary models not supported by existing spectrum calculators. 
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1 Introduction 

The era of the Large Hadron Collider is expected to be a data driven era. This statement is 
strengthened by a mix of numerous non-collider experiments, such as cold dark matter searches 
and flavor-centric probes, and near-future B factories. As experimental data becomes available in 
the form of measurements and null results, the value of the ability for rapid model discrimination 
increases. The SUSY Les Houches Accord (SLHA) is a protocol for interfacing the input and 
output of numerical tools that perform calculations within supersymmetric models |Tj. From 
here, all usage of the acronym "SLHA" will refer to the updated SLHA2 format and beyond. 
SuperLFV is a contribution to the growing library of SLHA tools, with an initial primary focus 
on calculating low-energy observables with charged lepton flavor- violation (LFV) [2^. 

The standard model (SM) with massless neutrinos prohibits LFV, due to an accidental sym- 
metry - accidental in the sense that no known reason protects the symmetry. Generically, ex- 
tensions of the SM will break this symmetry and exhibit new sources of flavor-changing neutral 
currents (FCNC). While it is observed that quarks exhibit FCNC and neutrinos oscillate flavor, 
thus far, charged leptons have not revealed flavor violation. The fact that neutrinos oscillate 
establishes charged LFV as a theoretical prediction. However, this prediction, given reasonable 
assumptions of unmeasured neutrino parameters, yields an undetectable /i — > 67 branching ratio 
of several orders beyond 10~^°. Hence any detection oi jj, ^ ej is an unmistakeable signature 
of new physics. For this reason, despite whether the branching ratio BR(^ — > ej) is measured 
or a stronger upper limit is placed, BR(/i — > 67) should be considered a standard benchmark for 
evaluating models. 

Supersymmetric models remain a primary focus of the field and, like any extension of the 
SM, generically exhibits LFV. The first calculation of the /i — >■ 67 branching ratio in a super- 
symmetric model was demonstrated in ref. ,3j. As will be emphasized (Section 3) LFV requires 
a source, i.e., any operator that couples lepton generations. For example, it is compelling that 
supersymmetric models with a neutrino seesaw mechanism, given reasonable parameters consis- 
tent with experiment, simultaneously provide an explanation of why /i — > 67 has not yet been 
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Table 1: Rates calculated by SuperLFV . Current experimental limits are listed at the 90% 
confidence level [7 . 



observed and a prediction that /i — > 67 may be observable with current technology. LFV decays 
for this particular model was first calculated in ref . [1] and comprehensively calculated in ref. |5] . 
The MEG experiment is a dedicated /x"*" — >■ e"'"7 experiment and has been acquiring data since 
2009 [6]. Thus far, the null result of MEG has provided a 90% confidence level upper limit on 
the — )• e+7 branching ratio of 5.7 x 10"^^. Consider the flavor violation to be manifested in 
a 1-loop diagram with a chargino and sneutrino with a mass-insertion (to|)i2. Then using an 
estimate from naive dimensional analysis 

i?fl(M^e7)-a^ ^^"^fif , (1) 

the MEG result implies that the off-diagonal entry (to? )i2 of the left-handed slepton mass matrix 
should not be larger than roughly 0.06% of the slepton masses. Here, ms is the typical sparticle 
mass, assumed to be 150 GeV to accommodate an MSSM interpretation of the muon anomalous 
magnetic moment measurement. This stringent constraint highlights the need for tools to rapidly 
calculate LFV observables. 

MEG's /i — > 67 search is just one of many near future and potential experiments, r — >■ 67 
and r ^ /i7 will be probed at the LHC and future _B-factories [5] d] . At minimum, these 
experiments should probe another order of magnitude in the various branching ratios; see Table 
^ Additionally, there are two other classes of LFV experiments — namely, the — >■ Se^ and 

N — >■ e~ N processes. The MUSIC project is a high-intensity muon source currently being 
constructed at Osaka University (Japan) and has proposed to extend the ^ — >■ 3e reach by 3 
orders of magnitude jTUj. Also, the PSI (Switzerland) may upgrade its existing muon source that 
provides 10^ muons/sec to 10^ muons/sec. If so, the proposed /i3e experiment may improve the 
/i — > 3e reach by 4 orders The tau variants, r — > 3e and t — !■ 3/x will also be further probed 
by the LHC and future i3-factories. 

Another class of LFV experiments is the search for neutrinoless decays of ground state muonic 
atoms, referred to as /i — e conversion in atomic nuclei ( "muon conversion," for brevity) . By 
expectations from the SM, when the muon decays, either a W~ boson is emitted outwards 
producing ordinary three-body muon decay or the W~ is captured by the nucleus, lowering 
the proton count by one ("muon capture"). Muon conversion seeks neutrinoless muon decay, 
with a free electron emitted. Coherent muon conversion refers to interactions that involve all 
nucleons, hence the amplitude scales with the number of nucleons. Coherent muon conversion 
requires an unexcited nucleus in the final state. There are four potential future muon conversion 
experiments, DeeMe, COMET, Mu2e, and PRISM/PRIME [11 [THIl [T^ . Each target atom 
provides a unique observable, as interactions vary for different nuclei. 

SuperLFV calculates the amplitudes of observables listed in Table[l]at one-loop level. SuperLFV 
accepts an SLHA spectrum file, which contains the couplings and mass parameters of the MSSM 
sparticles, and outputs the LFV rates. While this functionality is also included in related tools, 
SPhenodS], SuSeFLAV [T7], and partially by SUSY_FLAVDR iS], those tools currently do not accept 
externally generated SLHA spectrum files. SuperLFV is intended to be universal in the sense 
that it may be used with any tool that generates a SLHA spectrum file. 
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This separation of the spectrum generator and the observables calculator has the following 
merits. The source of LFV is allowed to be arbitrary. To date, all available SLHA-compliant 
spectrum generators use the supersymmetric seesaw mechanism as the source of LFV. However, 
grand unified theories [12] and flavor violating I?-terms from decoupled Z' models |i20j are ex- 
amples of other LFV sources that do not introduce new low scale particles. The SLHA spectrum 
file that SuperLFV uses as inputs may be generated by code or hand. Also, this approach facili- 
tates redundant checks of existing tools and allows focus on LFV-specific future enhancements to 
SuperLFV . Lastly, using an SLHA spectrum calculator, the SLHA allows one to specify arbitrary 
boundary conditions at a high scale (Planck, GUT, etc.) for an arbitrary model. The spectrum 
calculator will then evolve all parameters to the low scale, where SuperLFV may compute the 
LFV observables. 

Another feature of SuperLFV is the option to report all rates in terms of classes of contributing 
diagrams - i.e., contributions to 67 from diagrams involving charginos vs. neutralinos may 
be reported separately, as well as their interference contribution. Furthermore, any relevant 
coupling constants, effective and tree-level, are also optionally reported with a similar accounting 
for various contributions. Hence, SuperLFV is also designed to provide insight to the analytical 
behavior of the dynamics within the various LFV rates. 

This article is organized as follows. Relevant conventions are defined in Section 2. The 
physics of the included observables and their explicit expressions are reproduced in Section 3. 
In Section 4, the internal calculations of the code are detailed. All approximations used are 
disclosed in Sections 3 and 4 and Appendix A. The installation and operation of SuperLFV is 
documented in Section 5. 

2 Lagrangian and conventions 

All sign conventions, normalizations, conjugations, and matrix definitions are made to match the 
SLHA conventions, though the actual notation used may differ superficially in symbol choice. 
Irrespective of the high scale sources that generate the input parameters, SuperLFV assumes 
only the MSSM fields, including three Majorana left-handed neutrinos, as its complete effective 
theory. For clarity, the relevant notation used throughout this article is defined here. 
The superpotcntial is written 



where £12 = -1-1. Superfield notation is suppressed. Per the SLHA, the quark and lepton 
multiplets are written in the super-CKM and super-PMNS bases, in which quarks, charged 
leptons, and neutrinos are written as their mass states. That is, the sfermion components of the 
superfields have undergone the same rotations as their fermion counterparts. The Higgs vevs 
are defined as = w„/%/2 and (iJ^) = such that v'i + v'l^ (246 GeV)^ Accordingly, 

Yukawa couplings are defined at tree-level with \f2 factors; ^■g.^ = y^w^/v^. The gauge 
couplings of the hypercharge, weak, and strong interactions are denoted gy , g2i and (73. In all that 
follows, it is assumed that all couplings and mass parameters are evaluated at the supersymmetry 
breaking scale ("soft scale," for brevity); e.g., rusoh = ^/ttHJ^K^ or any arbitrarily defined scale. 

The SLHA requires that all running couplings and masses are defined in the DR scheme. 
There are some exceptions to this rule in defining the SMINPUTS block. However, the non-DR 
parameters of SMINPUTS are only used by spectrum calculators. As an observables calculator, 
only the output of those spectrum calculators are relevant to Super LFV. e.g., A spectrum 
calculator requires a'^^{'mz) defined via MS in the SMINPUTS block, but the spectrum calculator 
is required to output gauge couplings in the GAUGE block in the DR scheme, which is then used 
by SuperLFV . Hence, all SuperLFV running parameters are defined via DR. 

The soft parameters arc written 



W EE 



(2) 
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+ {-HuQauUR + HdQadd*ji + iJdZaeCfl + h.c.) 
+ ml^HlHu + m%mHd + {bH^Hd + h.c.) 



(3) 



All fermions are written in 2-component Wcyl notation, unless specified. Flavor indices are sup 



pressed, though symbols are arranged in matrix multiplication order; i. e., Q^m^ 



Isospin indices have been suppressed. The trilinear couplings follow the same isospin conven- 



Higgs soft mass operators 
H*^H-, H*Hd = H*^H2, and 



tions as their Yukawa counterparts; i.e., Hy^QauU*j^ = H^^afiQ 
have suppressed the following suppressed indices: H*Hu = 
HuHd = eafiH^H^. 

For the calculation of observables, all couplings are understood to be in the mass basis, 

which require mixing parameters. The remainder of this section will define the conventions used 
in determining masses and mixings. Starting with neutralinos, the mass operators are written 
as 





Ml 





-gyVd/'i. 






B 







M2 
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-9YVd/'2. 











^d 
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H° 



h.c. (4) 



Gaugino spinors have been redefined by the convention —iA — > A, siich that the off-diagonal 
2x2 sub-matrices are real. The unitary transformation matrix N that yields the neutralino 
mass states is defined by = Nab^b where A = [B, W^, H^, H^]. Throughout this article, 
capital indices A and B denote neutralino or chargino indices. By the SLHA, all sparticle mass 
indices are ordered by absolute value, with index 1 representing the lightest state. The SLHA 
allows for signed (real) masses for neutralinos and charginos. 
Chargino mass operators are written 



M2 g2Vu/V2 

giVdl\f2 M 



HI 



h.c. 



(5) 



The unitary transformation matrices U and V that yield the chargino mass states xf are defined 
by xT = UijXj and xt = 

To obtain the masses and mixings of sfermions in their mass states, it is conventional to form 
6x6 mass matrices for charged sleptons, up squarks, and down squarks and a 3 x 3 mass matrix 
for left-handed sneutrinos; i.e., 



UL 
Ur 



djL 
dR 



+ [44, 



]M? 



i>lMli>L. (6) 



Flavor indices have been suppressed, though ul = [uL,CL,iL], for example. Each 6x6 mass 
matrix of sfermion type / has the form 



2 2t 
'^)lR "^fRR 



(7) 



For squarks, 



uRR 



m 



uLR 
2 



dLL 



2 t I |2 



%/2 



w| + -\ydVd\^ + Dd^, 



(8) 
(9) 
(10) 

(11) 
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^Lk = + \\vM^ + D^,, (12) 

^Ilr = ^i.<^dVd-ydVuiJ.*)- (13) 



Again, flavor indices have been suppressed. For Yukawa matrices, = y\yu, despite being 

defined to be real and diagonal. The D-term for sfermion type / is 

Df = \{g^yY-glv'){vl-vl)l (14) 

where Y is the hypercharge generator and is the third SU(2) generator. Both generators are 
to be replaced by the corresponding eigenvalues of the fields they act on. 1 is the 3x3 identity 
matrix. Hypercharge assignments are defined via q = + Y, where q is th(^ electric charge of 
the representation of the sfermion. That is, q for left-handed representations of right-handed 
multiplets are opposite the charge of its Direct fermion component. The generator is defined 
as half the Pauli spin matrix a^. 



For charged sleptons, 



TOgflfl = rnf + ^\yeVd\'^ + Der,, (16) 
1 



1 
1, 



"^i^,z, = ml + -IveVdl"^ + De^, (15) 



^eLR = -^{ctlvd-vlvuH*). (17) 



Except for sneutrinos, sfermions of type / are rotated via a 6 x 6 matrix Rr, i.e. 



fi = {Rf)iJ 

where 



1^ 
Jr 



(18) 



RfMfRi (19) 



yields a diagonal matrix with the mass eigenvalues as the diagonal entries. 
The sneutrino mass matrix is written 

Mi = Ul^^s^lUpMNS + D,^ . (20) 

This expression assumes that m| is in the supcr-PMNS basis in wliic;h both fcrmionic and 
scalar components of the charged lepton supcrficld arc both rotated together, such that charged 
leptons are in their mass states, i.e., e^ V^e^ and U^eR with the same done to the 

scalar components. In the super-PMNS basis, the left-handed sneutrino states are rotated using 
the same transformation that rotates left-handed neutrinos to their mass states, i.e., vl Vi^vl 
and i>L — > V^Dl- Through the weak interaction operator W^leLi'y^vLi], the PMNS matrix is 
defined as {/pmns = K^K- 



3 Branching ratios 

In this section, the physics and expressions for the included LFV rates are reviewed. For each 
process, expectations for analytical behavior is discussed, as well as any approximations or 
omitted amplitudes in the SuperLFV code. 

LFV requires a source. Naturally, the LFV source will control the analytical behavior of LFV 
observables. Mainstream supersymmetry breaking models are generally designed to be free of 
observable charged LFV. This is due to the flavor problem of supersymmetry and the assumed 
irrelevance of flavor violation in probing the origin of electroweak symmetry breaking. The source 
is typically heavy particles, though light states that, say, couple only to third generation leptons 
are also conceivable. 
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We now commit to the MSSM in aU that follows. High scale dynamics may generate ofF- 
diagonal soft masses in the left-handed slepton mass matrix m?, the right-handed slepton mass 
matrix m^, or the matrix of trilinear couplings Ge- This allows the mass states of leptons 
and sleptons to exhibit LFV at tree-level interactions with neutralinos and charginos. Below 
the soft scale, the primitive LFV couplings will generate effective operators that only involve 
SM particles. The primitive couplings and any effective operators they generate are listed in 
Appendix A. 

The mainstream sources of LFV are neutrino seesaw models and grand unified theories 
(GUTs). The right-handed neutrinos of seesaw models will radiatively generate off-diagonal 
entries to the left-handed slepton mass matrix m| , due to Yukawa interactions with left-handed 
leptons. In this scenario, contributions to off-diagonal entries of the right-handed slepton mass 
matrix m| will be loop-suppressed. 

For SU(5) GUT models (with no seesaw mechanism), quark-mixing for the unified quark- 
lepton multiplets necessarily means lepton-mixing, induced by RG running between the Planck 
scale and GUT scale. When the SU(5) gauge group is spontaneously broken at the GUT scale, 
the MSSM as an effective model will be generated with off-diagonal entries to m|. Hence, 
7Ti| will be the dominant LFV source. These examples demonstrate that the LFV source is a 
model-building choice. While neither class of mainstream models yield the trilinear couplings 
{HdLaee\ + h.c.) as the dominant LFV source, it remains a possibility. 

In the following calculations, we follow the comprehensive work of ref. [5] with corrections to 
the Ci — )■ 3ej decay provided by ref. [21] . 

3.1 Cj — )■ ej7 

The effective lagrangian for ^jl* ^ defined at the scale of the lepton mass mg., is 

+ ^e,a,,F'^''{A^'^PL + A'i^PR)e, + h.c. (21) 

The electromagnetic form factors A^-^^ and A^-^^^ are defined to exclude the tree- level value —J'-', 
and do not contribute to on-shell — > 6^7, as made explicit by factoring out the photon's 
external momentum squared from the couplings. Au^ and Am will be relevant later to 
Ci — > SCj and muon conversion, where off-shell photons contribute. 

The Ci ej7 branching ratio definition deviates from the conventional definition of a branch- 
ing ratio in which a partial width is compared to the total width. Comparison of the par- 
tial width r(ei — ej7) to T{ei — > ejVjVi) is the standard observable, evaluated in the limit 
771^. <^ mg. <^ TOn/. For r decays, this definition deviates significantly from the conventional 
definition. To lowest order, the partial width of — >■ ejVjVi is 

r(e, ^ e^VjVi) = ^^^J . (22) 



Consequently, the effective theory given by eqn. (21 1 yields 



BR(e. ^ e,7) - 48^^^ (l^il' + • (23) 

The following remarks are regarding the expected analytical behavior of — >■ 6^7. 

• At I-loop order, the couplings A^^l ^^"^ generated by loops involving either 

charginos or neutralinos. These form factors generally scale linearly tan /3. This can 
be understood by considering diagrams involving all sparticles in their gauge-interaction 
states ("gauge states," for brevity). Then diagrams can be distinguished as those with 
only gauginos in the loops vs. those with higgsino-gaugino mixing in the loops. Then the 
higgsino-gaugino mixing manifests as mass-insertion involving the up- type Higgs vev Vu- 
Also, the Yukawa coupling of a lepton and down-Higgsino will contain 1/vd- Together, 
the tan/3 dependence is formed. 
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Experimental constraints require ofT-diagonal soft masses of charged sleptons to be small. 



Yet, ofF-diagonal masses for sneutrinos, in the Super-PMNS basis of eqn. (20), may be 
large. Therefore, diagrams involving sneutrino mixings should dominate over those with 
slepton mixing. Alternatively, chargino diagrams should generally dominate over neutralino 
diagrams. 



The tensor operators of eqn. (21) force a helicity change. Amplitudes that change the 
lepton helicity from right to left will involve the Yukawa coupling of the decaying (heavier) 
lepton, for diagrams with higgsino-gaugino mixing. This will dominate over diagrams that 
involve Yukawa coupling of the lighter final-state lepton. If the LFV source is m|, the 
initial lepton must be right-handed for the diagram to be proportional to larger Yukawa 
coupling. Hence, efn — >■ cljJ should dominate over CLi — >■ e^j^y when m| is the LFV 
source. Equivalently, the A^j^ coupling should dominate over A2j^. These statements are 
reversed if the LFV source is m|. 

Furthermore, if m| is the dominant LFV source, there are no charged winos in the loop 
diagrams. Only amplitudes with loops involving either a pure bino or a bino-Higgsino 
mixing remain. 

At one-loop, a non-holomorphic coupling of leptons to the up-type Higgs iJ^ generates 
the effective coupling of the form emyij Zij H'^* clj , if an LFV source exist for m^, m^, or 
o-e [22j . Here, Zij parameterizes the radiatively-generated coupling. This up- type Higgs 
couplings yields two more mechanisms for mediating LFV. Presently, neither have been 
included in SuperLFV . 

— Using this effective coupling at one-loop, the neutral Higgs bosons H^, and A" 
will mediate — 6^-7 |23| . However, such diagrams are doubly Yukawa suppressed. 
At two-loops, diagrams classified as Barr-Zee diagrams avoid this suppression and 
dominate. It can be demonstrated that for reasonably light sparticles, this Higgs- 
mediated LFV is negligible. It becomes dominant when the ratio of the sparticle 
mass scale to the mass of the pseudo-scalar Higgs msoft/^Tiyi grows beyond roughly 
40-50. 

— At tree-level, the Higgs propagator from the non-holomorphic coupling vertex may be 
connected to another lepton or nucleon to yield — 3ej 22J and muon conversion 
|24| . This is irrelevant to — > 6^7, but pointed out to avoid confusion, as both types 
of Higgs mediated LFV enter in other LFV observables. 

If trilinear couplings are the dominant LFV source, the dominant loop diagram (cast in 
mass states) will involve only a neutralino propagator and slepton propagator. A chargino 
would create a left-handed sneutrino. The trilinear coupling would flip the neutrino helicity 
to the decoupled right-handed sneutrino. Hence, charginos do not contribute to this pro- 
cess. The neutralino-slepton loop will be much like those with m? or m- mass-insertions. 
In the gauge-interaction basis, there will be two one-loop amplitudes: one with a pure bino 
and another with a bino-higgsino mixing. Either can be made to dominate, resulting in 
non-trivial behavior. 

— In a loop with a pure bino for the neutralino propagator, there will be a single mass 
insertion by the trilinear coupling of ivd{ae)ij/V2 along the slepton propagator. The 
tan/3 dependence of this amplitude is controlled by the Vd vev, or v/\/l + tan^ f3 . 
The left-right mixing of the trilinear coupling is sufhcicnt for the required helcity flip 
of the electromagnetic dipole operator. 

— In a loop with a bino-higgsino transition along the neutralino propagator, since the 
higgsino will cause one helicity flip, two left-right mass insertions are required along 
the slepton propagator. One trilinear should be the LFV mass-insertion ivd{cLe)ij 
The other mass-insertion should be the usual flavor-preserving left-right slepton mix- 
ing. Hence, these three helicity flips satisfy the electromagnetic dipole operator. 
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Each left-right mixing introduces one factor of the Higgs vev Vd to the amplitude. 
The bino-higgsino transition inlcudes a w„ factor, while the lepton-slepton-higgsino 
Yukawa coupling includes a 1/vd factor. This results in an overall VuVd factor, or 
v'^ tan/3/(l + tan^ (3). That is, this amplitude decreases with tan /3. Also note that, 
if the flavor-preserving trilncar couplings are negligible, this amplitude itself becomes 
neglibible. 

Unlike the m| and m| LFV sources, a LFV trilinear does not force the initial lepton to be 
left- or right-handed. Therefore, the amplitudes for em — > e^j^ and e^i — > enjj (or the 
magnitudes for A^j^ and A^^l) ^^'^ expected to be comparable for trilinear LFV sources. 

• Loops with purely SM particles, i.e. W bosons and neutrinos, are omitted in the SuperLFV code, 
as it is for all other rates to follow. In principle, neutrino oscillations will yield non-zero 
LFV. However, the smallness of the neutrino masses place contributions to the /i — > 67 
branching ratio on the order 0{m^/my[r)^j several orders below 10^^° for reasonable neu- 
trino parameters. 



3.2 ei^ejejel 

The amplitude for — >■ 3ej decay includes diagrams with photon exchange, Z exchange, Higgs 
exchange, and neutralino-slepton (chargino-sneutrino) boxes. Defined at the scale of the lepton 
mass TTT-e. , the effective lagrangian for — > Sej decay that provides photon exchange is 

2 

-CD ^e^[e,J,iA[iPL+A[^J,PM[e,re,] + ^-^[e^a,,q^ {A'i^PL + A^^Pn)e,][e,re,]. (24) 

Z exchange is governed by 

■ 2 

- £ D ^[e,7^(F^Pi + Fi^Pn)e,][e,r(Ze,PL + Z,,Pn)e,] (25) 

where is the tree- level Z coupling to a chiral fermion ?/;, given by 

Z-4,^tl-q^siT?e^. (26) 

Neutralino-slepton (chargino-sneutrino) box diagrams are governed by 

- £ D ie^B\'^ [eJJ^PLe^] [eji^'PLCj] + le^ B\'^[ej-f^PRe^] [e,7^Pfle,] 
+ie^B^^j^[ej-i^,PLei][ej-f^PRej] + ie^ B^^j^[ej'-if_,PRei][ej-i^ PlCj] 
+ie^ Bl^^[e.jPLe.j\[ejPLej] + ie^ B'-^j^[ejPRei][e.jPRej] 

+ie^Bll[eja^,PLe,][e,a''''PLe,] + le^ Bl^^[eja^,PRe,][e,a>''' Pne,]. (27) 

All form factors are supplied in Appendix A. 

As with e; — ^ 6^7, the — 3ej branching ratio is defined as the ratio of its partial width to 
the ei — > ejVjVi. This results in a branching ratio of [51 121) 



BR(er ^ e7e-e+) = :i^[|Ali|2 -f |Al^|2 _ 4Re(#iA^^^* + A^i#/^) 

Y^{\B-L? + \B':^H?) + m'iL? + \BU) 
Re{BliBl^*+Bl^^Bl^*) 

Ir^{A^,B^* + A^rB^; + A^,B^* + A^rB^*) 
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4 .... .... .... .... 

T? ( A^^ R^.?* _L A^J R*.?* _L A^J R^.?* _L A^J pU*" 



where 



1 

3'' 

4Re(A;i^^ + + 2Re(^;if^ + J^^F]^^ 

8Re(A^i^+A^i^)-4Re(^^J^+^^^^)] (28) 



F ■' 7 

P^L - . L ^\ . (29) 

Sm OwTTlf^ 

F^R = FI'lIl^R^ (30) 



n^R - . 2, \ : (31) 

sm Owtn^ 

FlL = ^^L^kofi. (32) 
The following is noted about expectations for the — > Se^ branching ratio. 

• On-shell photon exchange dominates due to the tan^ j3 enhancement in the branching ra- 
tio and a relatively large logarithmic enhancement (^In^^ ~ relative to all other 

couplings, which are all defined to be of dimension GeV^^ (for comparison with the Fermi 
constant Gp)- It can be verified with the SuperLFV code that the photon-exchange dia- 
grams dominates these decays, with all other couplings typically being two to several orders 
of magnitude smaller. Even in the case of low tan/3, the logarithmic enhancement of the 
dipole exchange contribution is sufficient to ignore all other contributions to percentile 
accuracy. This yields a quite rigid prediction of 

BR(ei ^ 3e,) _ 2a ( _ 11\ ^^^^ 



BR(ei — >■ 6^7) Btt \ TTie^ 8 
This provides numerical values of roughly 0.61% — >■ 3e), 0.22% (r 3/i), 1.1% (r 3e). 

• The Z-cxchangc diagrams can be classified as eight types, four each for neutralino and 
chargino loops. Within this division, two diagrams involve loops with a Z attached to 
either the neutralino (chargino) or slepton (sneutrino). The other two diagrams are LFV 
self-energy loops on either the initial lepton or the final lepton. The analytical behavior of 
these Z-exchange diagrams was elucidated by ref. [35] . These authors divide contributions 
into those with couplings Fji and (see Appendix A). Fji are all Yukawa suppressed 
diagrams. Next, they demonstrate that, in the absence of chargino-mixing, all contributions 
by a pure charged wino cancel exactly within Fi^. 

• As mentioned earlier, the non-holomorphic Higgs coupling is omitted in the photon ex- 
change and Higgs exchange amplitudes for — )• 3ej . Despite a tan^ /3 scaling, Higgs- 
exchange is unable to overcome Yukawa suppression for both muon and tau decays. In 
the most complete study of ei 3ej, it is demonstrated that, in the then-experimentally 
allowed parameter space of minimal supergravity (mSUGRA), Higgs exchange remains a 
few orders of magnitude smaller than photon exchange for tan /3 values as high as 50 [H] . 
Yet, under more extreme conditions of heavy sparticles, a light pseudo-scalar neutral Higgs 

, and large tan/3, Higgs exchange may compete with photon exchange. 

3.3 Muon conversion fi^ N N 

This process is governed by the same classes of diagrams as — )• Sej. For QCD concerns, 
only quarks lighter than the QCD scale are included the effective theory. Furthermore, strange 
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quarks do not contribute to coherent muon conversion for the foUowing reasons. In the limit 
of isospin symmetry amongst light quarks, there is no vectorial coupling to strange quarks in 
the nucleus. The axial-vector couplings to strange quarks yields an incoherent interaction due 
to spin transitions and is therefore neglected. At the quark level, the lagrangian for photon 
exchange is 

Q 

■ 2 

ip'^m _ ^ — - 

^-[e.a^Xi^'^LPL + A'iRPR)e^] ■ V q^^q] (34) 
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where the sum is over up and down quarks with electric charge g,. Z exchange is described 

by 

1.9| 



^ ^ 2 S^'"^'^^^^'^^ + ^« ^^^^^ + ^9Jb-7^9]. (35) 

Z q 

The neutralino-slepton (chargino-sneutrino) box diagrams are described by 



£ D ^[67^(^9.^^ + Dq,Pn)mrq]- (36) 

9 

The four-point effective couplings F^, F^, Dg^^ , and Dg^^ are supplied in Appendix A. Note that 
Higgs exchange diagrams are currently omitted. To calculate muon conversion amplitudes, the 
quark-level effective lagrangian is then used to create an effective theory at the nucleon-level 

m- 

The muon conversion branching ratio is defined as the ratio of muon conversion partial width 
to the muon capture partial width, where the final state nucleus of muon capture may be excited. 
That is, 

BRi,-N^e-N)^£,^^^-^^. (37) 

The coherent muon conversion decay width for a muonic atom of atomic number Z and N 
nucleons is [S] 

Tifi-N^e-N) = Aa''^\FM?ml[\Z{All-Al\)-{2Z + N)D.^,-{Z + 2N)DaA^ 

+ \Z{Al\ - All) - (2^ + N)D,^, - (Z + 27V)^,J2] (gg) 
where the Z-exchange and box contributions have been combined as 



^QL - ^qL + r,-2n 9 ' (.3yj 

z sm u^m-^ 



Dg, = Dg,\L^R. (40) 

Zes and |F/v| describe an effective nuclear charge and nuclear matrix element, respectively, and 
are defined in the pioneering framework provided by Weinberg and Feinberg [26] . Table [2] lists 
the values used in the code. 

To address the nuclear aspects, a few assumptions are made. Currently, SuperLFV calculates 
muon conversion using the general approach of Weinberg and Feinberg, applied to the MSSM 
by Hisano, et al. These assumptions are listed. 

1. For Z-exchange and box diagrams, only vector currents are expected to be relevant, since 
axial-vector quark couplings yield spin-dependent nuclear effects. This approximation is 
valid for sufficiently heavy isotopes. 

2. The densities of protons and neutrons are assumed to be equal and constant. 

3. The muon is treated non-relativistically, and the muon wavefunction overlap with the 
nucleus is taken to be a constant. 
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gpff \Fn\ r{pN^Uf,N') R 

ilAl 11.62 0.64 0.7054 X lOS"^ 351 
tlTi 17.61 0.535 2.59 x 10*^8-^ 203 

Table 2: Nuclear data for muon conversion. Here, Z^s and \Fi\f \ are the effective nuclear 
charge and a nuclear matrix element that describes charge distribution, respectively, of eqn. (38). 
Estimates for Zcff and \Fn\ are from |28| . Estimates for the experimental muon capture partial 
widths T{fiN -> i/^iV') are from [23 



4. _FV is a weighted overlap of the muon wavefunction and the nucleons, and hence applies 
the approximations of the previous two points. 

5. The final state electron is treated non-relativistically as a plane wave, undistorted by the 
nuclear electric charge. 

6. The electron mass is neglected. 

The state-of-the-art technique of addressing the nuclear aspects is currently authored by 
Kitano, et al. Those authors apply the approach of Czarnecki, et al. (3D], which includes 

an improved scheme for nuclear considerations of off-shell photon exchange. These groups also 
solved for the initial-state muon and final-state electron wavefunctions via the Dirac equation. 
Kitano, et al. also applied various schemes for nucleon distributions and experimental neutron 
distributions obtained from pionic atoms. They find that coherent muon conversion rates for 
heavy elements can differ by ~ 20 — 30% from the method of ^\ . 

The following points are noted about expectations for the muon conversion rate. 

• For much of the MSSM parameter space, the photon-exchange amplitudes dominate the 
muon conversion rate. In this case, there is roughly a fixed correlation of order a between 
a muon conversion rate and BR(// — > 67) that only depends on the atomic element. For 
example, using the framework of Weinberg and Feinberg and the nuclear approximation 
scheme of Bernabeu, et al., one arrives at a relation 



BR(^- N ^ e- N) 



where values for R are listed in Table [2j This convenient guideline breaks down in the 
limits of both small and large tan /3, as noted below. 

Unlike — ?► 3ej, on-shell photon exchange is not boosted by a large logarithm. For small 
tan P and ^ < 0, the on-shell photon exchange contribution may have large cancellations 
amongst its gauge state amplitudes [5]. 

• Also unlike — ?> 3ej, Higgs exchange is not doubly Yukawa suppressed. The neutral 
Higgs bosons couple proportional to the nucleon mass, via the superconformal anomaly. 
For this reason and the previous point, Higgs exchange may be demonstrated to dominate 
for experimentally accessible parameter space. For example, for mSUGRA parameters of 
unified scalar and gaugino soft masses of 1 TeV at the unification scale, i.e. mo = mi/2 = 
1 TeV, muon conversion for aluminum begins to deviate from photon exchange dominance 
for a heavy Higgs lighter than 400 GeV [21]. 

4 Numerical calculations for the LFV observables 

In this section, the procedure of the numerical calculations is made transparent. 
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4.1 Calculation procedure 



An SLHA spectrum file supplies all input for the calculation of amplitudes. For branching ratios 
evaluated at low energy, couplings and masses that arise due to solely to kinematics are evaluated 
using hard-coded values compiled by the Particle Data Group (PDG); e.g., a(0) or pole masses. 
The values of all Yukawa and gauge couplings and running mass parameters are read in at the 
scale supplied by the SLHA "Q=" convention. For example, immediately following the header 
BLOCK GAUGE Q= 1 .OOOOOOOOE+03, gauge coupHngs are supplied at the RG scale of 1 TeV. 
These running parameters are not loop-corrected by SuperLFV , since they may have already 
been loop-corrected by a spectrum calculator. That is, SuperLFV accepts and uses parameters 
as is. 

The Higgs fields do not enter in any of the calculations currently offered by SuperLFV . 
Therefore, their masses and mixing parameters are ignored. 

Next, the sfermion mass matrices, along with all other required inputs [e.g., those necessary 
for left-right mixing and D-terms) are supplied by the SLHA spectrum file in the DR scheme and 
used to calculate the running sfermion masses and mixing angles at the (3= scale. Parameters 
that form the neutralino and chargino mass matrices are read in. The physical masses and 
mixings are then calculated. Wavefunction renormalization of fields is currently not performed 
in any diagonalization process. 

Naively, it would seem simpler to use the sparticle mixing matrices supplied by an SLHA 
spectrum file. However, as allowed by the SLHA, those mixings may not be calculated or 
reported in same manner as the mass matrices. For example, the masses may be radiatively 
corrected while the mixing matrices are not. This matter is very important to calculations 
involving snuetrinos and charged sleptons, as they tend to include nearly degenerate states. 
The orthogonality and near-degeneracy of physical sneutrino and slepton states cause large 
cancellations within the photon penguin diagrams, hence the parameters used must be sufficiently 
precise. Also, insufficient accuracy by an SLHA spcictrum calculate may interchange rows in the 
mixing matrices, leading to garbage results. Since SuperLFV is not tied a specific spectrum 
calculator, this uncertainty prohibits usage of mixing parameters supplied by SLHA spectrum 
calculators. To ensure the integrity of calculations, SuperLFV calculates its own physical masses 
and mixings. 

The SLHA allows for multiple entries of the same BLOCK at different scales. For exam- 
ple, it is common that the BLOCK GAUGE is supplied at the GUT scale and soft mass scale. If 
SuperLFV encounters a block name that has already been read in, the existing block is overwrit- 
ten; i.e., the block further down the file is retained. The code checks that the defining scales 
of Yukawa couplings, gauge couplings, and mass matrices all match. If they do not match, the 
code aborts. 

Next, using the running masses and mixings evaluated at the soft scale, primitive couplings 

of mass states (such as snuetrino-lcpton-chargino couplings) are evaluated. Then, also at the 
soft scale, all relevant effective couplings are then evaluated using only the running masses and 
primitive couplings. 

Finally, observables are calculated at the soft scale and supplied as output. In principle, one 
should run each effective coupling c down to the scale of the initial state lepton mass Wj using a 
full RG evolution or a leading log approximation, Ac = /3cln(TOi/TOsoft)- Since the electromag- 
netic tensor couplings A^^l ^^'-'^ ^2R tend to dominate LFV processes in supersymmetric models, 
SuperLFV will RG evolve them from the scale Wsoft at which they are generated to the muon 
mass m^. This RG running is performed via the leading-log approximation. That is. 



where A represents one of the A^^ and A]^^^ couplings. In our context, other corrections to other 
operators are small and neglected. 




(42) 
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SuperLFV 



SPheno 



LFV source BR(/i — > 67) BR(/i — > 3e) ratio BR(/i — >• 67) BR(/i — > 3e) ratio 

(m|)i2 = (lOOGeV)^ 4.74x10"^^ 3.27x10"^'' 0.69% 3.57x10"^^ 2.19x10"^"* 0.61% 
(mi)i2 = (lOOGeV)^ 5.24 x 10^^^ 3.47 x 10^^^ 0.66% 1.53x10"" 8.84x10"^^ 0.58% 



Table 3: Numerical comparison of SuperLFV vl.O and SPheno v3.2.1. For illustrative 
purposes, an mSUGRA model augmented by an LFV source is presented. The mSUGRA model 
is parameterized as niQ = 2.5 TeV, mi/2 = 1.5 TeV, Aq — 0, tan/3 = 10, and /i > 0. Here, 
"ratio" is defined as BR(/i 3e)/BR(/i ej). 

4.2 Comparison to existing SLHA tools 

The current version of SPheno readily outputs an SLHA spectrum that meets the minimum re- 
quirements of SuperLFV [^Furthermore, SPheno currently calculates the branching ratio BR(/i — > 
67), along with other variants. This allows for convenient comparison. It is worth repeating that 
SPheno 's ability to calculate LFV observables is not fully redundant with SuperLFV . A major 
impetus for SuperLFV is to divorce SLHA spectrum calculators from the calculation of LFV 
observables, allowing for more robust model selection. 

Procedural comparisons will be made first. The electromagnetic dipole couplings 
A^2Ri when squared for observables, significantly lower the observable. For ^ — 67 and mgoft 
between 100 and 1,000 GeV, the rate BR(/^ 67) is lowered by about 12% to 17% [31]. SPheno 
currently does not perform this correction. Also, SPheno also performs a loop-corrected wave- 
function renormalization of the gaugino and higgsino fields; SuperLFV does not. Neither tool 
performs a proper matrix treatemcnt for the wavefunction renormalization of the squarks and 
sleptons. 

Numerical comparison examples for SuperLFV vl.O and SPheno v3.2.1 are listed in Table [Sj 
All discrepancies have been traced back to the neutralino and chargino mixing matrices. All 
running couplings, running masses, and sfermion mixing for SuperLFV and SPheno have been 
carefully checked to match to several digits of accuracy, typically 1 part in 10^ at minimum. 
The fact that SPheno performs a wavefunction renormalization of the chargino and neutralino 
matrices while SuperLFV does not, creates a small discrepancies in the neutralino and chargino 
mixing matrices of SPheno and SuperLFV . The near-degeneracy of sneutrinos and charged slep- 
tons, along with the orthogonality of their states, results in precise cancellations in calculating the 
electromagnetic dipole couplings A^2l ^'^'^ -^in- Small deviations in the neutralino and chargino 
mixing matrices of SPheno and SuperLFV will result in imprecise cancellations. Given that rates 
are proportional to |^2lP + I ^2^1^' "^^^^ effect is that SuperLFV 's rates are higher than 
Spheno's but on the same order. Note that this is only an artifact of using a SPheno spectrum 
file as input for SuperLFV . 

Given a spectrum over which an author has full control, SuperLFV 's precision is well- 
controlled. To demonstrate this, the tan/3 enhancement of /i — 67 is demonstrated in Fig. [T] 
Recall how this enhancement arises. The Vu factor of tan/3 arises from the higgsino-gaugino 
mass-insertion within the /i — > 67 loop. The l/vd factor arises from the muon-higgsino Yukawa 
coupling. Vu is encoded into the neutralino (chargino) mixing matrix. Vd is encoded into the 
muon-neutralino (chargino) Yukawa couplings, which is composed of elements from the mixing 
matrices of the slepton (sneutrino) and neutralino (chargino). Hence, without proper accuracy, 
the tan/3 enhancement is easily destroyed. 

^In principle, other SLHA spectrum calculators may also meet SuperLFV 's minimal spectrum requirements, 
however most have omissions in the output - a trivial obstacle for the authors of SLHA spectrum calculators to 
resolve. 
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5 Installing and operating SuperLFV 
5.1 Installation 

SuperLFV is packaged as a single file, named superlf v-1 .0.0. zip. It is installed by the following 
process. 

1. Obtain the SuperLFV code from the HepForge distribution service at: 



http : // superlf V . hepf orge . org 



2. Double-click superlfv.zip. On most modern computing platforms, a .zip file will auto- 
matically expand into constituent files. A resulting directory named SuperLFV will appear. 

3. In a terminal window, change the directory to SuperLFV. 
5.2 Immediate usage 

SuperLFV is supplied as a single pre-compiled Java JAR file superlf v. jar which, aside from an 
SLHA spectrum file, contains all that is necessary to start using SuperLFV immediately. There 
is no need to expand the JAR file into constituent files. As Java code, SuperLFV should run 
on any modern system with a Java Virtual Machine. An example default input file is supplied 
as input. spc. To immediately run the code using the default input file, using a command line 
interface type: 

Java -classpath superlf v. jar superlf v 

This command is tedious and may be vastly simplified, as described below. By default, only the 
LFV rates are displayed. The supplied example input. spc yields the following output. 




Figure 1: A demonstration of the tan/? enhancement of — 67. All models (A, B, C) 
shown have spectra generated by mSUGRA boundary conditions. In common, they have ^0 = 
and /X > 0. Model A represents toq = 2 TeV, mi/2 = 1.6 TeV, and (m|)i2 = (100 GeV)^. Model 
B represents mg = 500 GeV, mi/2 = 300 GeV, and (rn|)i2 = (3GeV)^. Model C represents 
mo = 1 TeV, mi/2 = 800 GeV, and (m|)i2 = (lOGeV)^. 
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c Displays the primitive and effective couplings used for calculations. 

C Same as c but with the contributions of various amplitude classes shown, 

i Allows an input file to be specified. 

p Displays the parameters (running masses and mixings) used for calculations, 

s suppresses the output of observables. 

w Displays the observables with contributions from various amplitude classes. 
Table 4: A list of command line options. 



BR(mu- -> e- gamma) = 1 . 891993318E-14 
BR(tau- -> e- gamma) = 8 . 773022513E-36 
BR(tau- -> mu- gamma) = 1 . 925024531E-34 

Observable: R(mu N -> e N) 

BR(mu Ti -> e Ti) = 1 . 197657937E-16 
BR(mu Al -> e Al) = 7 . 041410720E-17 

Observable: BR(1 -> 1' 1' 1') 

BR(mu- -> e- e- e+) = 1 . 311991365E-16 
BR(tau- -> e- e- e+) = 1 . 734152699E-38 
BR(tau- -> mu- mu- mu+) = 9 . 690935025E-38 



5.3 Long-term usage 

Invoking SuperLFV is vastly simplified by creating an alias. On Unix variants, one should edit 
the ~/. profile file to include the following line. 

alias superlfv=" Java -classpath (iirector?/-pai/i/superlfv. jar superlfv" 

Here, directory-path^^ should be replaced by the directory path that contains the superlfv. jar 
file. This may be obtained by typing "pwd" in a terminal window of a Unix variant system. To 
invoke the change made to the ^7 .profile file, either open a new terminal window or manually 
invoke it via "source ~/ . profile" . From here on, SuperLFV may be invoked with the following 
simple command: 

superlfv 

All example commands to invoke SuperLFV will use this simplified command. 

5.4 Options 

Table [4] shows a list options that are supplied in the following format. 

superlfv -options [input file] 

These options may be used in conjunction. For example, to display the spectrum parameters 
used for the calculation in addition to supplying an input file different from input . spc, one may 
use: 

superlfv -pi myinput . spc 

The ordering of options is arbitrary. Options -ip would produce the same results in this example. 

One feature of SuperLFV is to report quantities broken down into various contributions. For 
example, using the option -w, the n 3e branching ratio will explicitly show contributions from 
photon-exchange, Z-exchange, box diagrams, and interference terms. The same insights may be 
gained from primitive and effective couplings using the -C option. 

Output using the -w apears in the manner below. 
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Observable: BR(1 -> 1' gamma) 

BR(mu- -> e- gamma) = 1 . 891993318E-14 = 4.863598412E-16 (neutralino) + 2 . 547321573E-14 (chargino) 

+ -7.039642383E-15 (interference) 
BR(tau- -> e- gamma) = 8 . 773022513E-36 = 8 . 773022513E-36 (neutralino) + 0.00000000000 (chargino) 

+ 0.00000000000 (interference) 
BRCtau- -> mu- gamma) = 1 . 925024531E-34 = 1 .925024531E-34 (neutralino) + 0.00000000000 (chargino) 

+ 0.00000000000 (interference) 



Output using the -C apears in the manner below. 



Couplings: lepton-photon-lepton 1-loop vector and tensor couplings 

A_1R-11 = 9.604168151E-09 = 1.283972049E-09 (neutralino) + 8.320196103E-09 (chargino) 
A_1R-12 = -6.564753992E-14 = -1 .765258224E-14 (neutralino) + -4.799495768E-14 (chargino) 
A_1R-13 = 8.132763928E-28 = 8. 132763928E-28 (neutralino) + 0.00000000000 (chargino) 
A_1R-21 = -6.564753992E-14 = -1 .765258224E-14 (neutralino) + -4.799495768E-14 (chargino) 
A_lR-22 = 9. 60492320 lE-09 = 1.284201890E-09 (neutralino) + 8.320721312E-09 (chargino) 
A_lR-23 = -2.353389164E-24 = -2.353389164E-24 (neutralino) + 0.00000000000 (chargino) 



Output using the -p apears in the manner below. 



SuperLFV 1.0 

Scale q = 1000.0 

Higgs parameters : 

tan (beta) = 10.00000000000 
v_u = 242.39690736872 
v_d = 24.23969073687 
V = 243.60587700000 
mu = 399.82391000000 

Gauge couplings: 

g_Y = 0.36283193300 
g_2 = 0.64585703200 
g_3 = 1.07837289000 

Yukawa couplings: 

y_u = 8.493623820E-06 

y_c = 0.00359564127 
y_t = 0.87077568600 
y_d = 1.364646190E-04 
y_s = 0.00286579735 
y_b = 0.13692728200 
y_e = 2.982727560E-05 
y_mu = 0.00616732139 
y.tau = 0.10376481700 

Fermion mixing matrices: 
V_CKM = 

0.97419203200 0.22569459200 0.00344018712 
-0.22564536100 0.97335925700 0.04069308400 
0.00583567101 -0.04041914040 0.99916577100 

U.PMNS = 

1.00000000000 7.686990450E-08 0.00000000000 
-7.686990450E-08 1.00000000000 0.00000000000 
0.00000000000 0.00000000000 1.00000000000 



Spectrum: Squarks 
Masses : 
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Mixing matrices: 
U.u = 

1.971336933E-05 2 . 332759829E-04 0.38641862210 -5 . 572092936E-07 -6.437877395E-06 0.92232347560 
-1.138264522E-04 -0.00134740573 -0.92232253487 2 . 739826674E-06 2 . 970168476E-05 0.38641857140 
-0.07314197920 -0.99728970343 0.00133727625 -1 . 372091198E-06 -0.00784958826 -3 . 065230137E-04 
0.99732153150 -0.07313980192 -1 . 483096383E-05 1 . 854684657E-05 -5 . 758580275E-04 3 . 391928336E-06 
-1.848175312E-07 0.00787063396 -4 . 037200977E-05 6 . 224076218E-06 -0.99996902521 7 . 943845986E-06 
-1.859720289E-05 -5.702575808E-08 2 . 744681143E-06 0.99999999980 6 . 223738233E-06 -5.453256277E-07 

U_d = 

-0.00564646735 0.03910914089 -0.99675023977 3 . 251631993E-06 2 . 938364194E-05 -0.07019651810 
0.00184956358 -0.01281420614 0.06973474530 -2 . 145209511E-06 -2 . 886612624E-04 -0.99748150720 
-1.397029538E-06 -0.00699365618 -3 . 289657427E-04 2 . 846389612E-05 -0.99997542612 3 . 562264875E-04 
3.336819281E-04 -6 . 002136498E-07 1 . 450037733E-06 0.99999994392 2.846733857E-05 -1 .431055699E-06 
-0.16598733285 0.98522716484 0.04031152804 5 . 610198104E-05 -0.00690715109 -0.01014434873 
-0.98610992821 -0.16608670521 -9 . 472630988E-04 3.289159902E-04 0.00116336550 2 . 386035123E-04 



Spectrum: Charged sleptons 
Masses : 

m2_el = 6.052707465E04 
m2_e2 = 6 . 376245232E04 
m2_e3 = 6 . 377409414E04 
m2_e4 = 4 . 160578137E05 
m2_e5 = 4.221118179E05 
m2_e6 = 4.221332878E05 



Mixing matrices: 
U_e = 

-9.110035956E-19 3 . 003274668E-15 0.99978452461 4 . 591381300E-17 1 . 865954857E-16 0.02075823592 
0.08621816137 -0.99627553490 4 . 956493417E-15 5 . 103829867E-07 -0.00121951244 1 . 038557432E-16 
0.99627628127 0.08621809678 -4 . 274129627E-16 5 . 897817520E-06 1 . 055405394E-04 -8 . 956081394E-18 
7.633941527E-22 4.531569693E-18 -0.02075823592 -1 . 251423047E-16 -2 . 967572004E-16 0.99978452461 
-3.415966935E-09 -0.00122407083 4.013172943E-17 -3.375193690E-10 0.99999925083 1.456935011E-16 
5.919859989E-06 1 . 693127777E-11 2. 154393278E-17 -0.99999999998 -3.374799962E-10 -2.541013616E-17 



Spectrum: Charginos 
Masses: 

m_Cl = 180.67092936012 
m_C2 = 420.58024290262 

Mixing matrices: 
U = 

0.93957277159 -0.34234924694 
0.34234924694 0.93957277159 



V = 

0.98331246394 
0.18192470493 



-0.18192470493 
0.98331246394 



Spectrum: Neutralinos 
Masses : 

m.Nl = 100.87739321550 
m_N2 = 181.06276934770 
m_N3 = -405.56812343355 
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m_N4 = 420.17616587035 

Mixing matrices: 
N = 

-0.98936250240 0.04879451751 
-0.08376027510 -0.96033420653 
0.05528440736 -0.08366717812 
-0.10530759453 0.26149017215 



-0.12960484767 0.04453669689 

0.23511433667 -0.12435303281 
-0.69738290128 -0.70964817100 
0.66451555373 -0.69206378265 



Spectrum: Sneutrinos 
Masses : 

m2_nul = 5.461111497E04 
m2_nu2 = 5 . 769697773E04 
m2_nu3 = 5 . 770810200E04 

Mixing matrices: 
U_nu = 

0.00000000000 0.00000000000 1.00000000000 
0.09026203273 -0.99591805157 0.00000000000 
0.99591805157 0.09026203273 0.00000000000 



Trilinears : 
a_u = 

-0.00499689032 2 . 262252780E-09 3 . 249917170E-08 
9.576890890E-07 -2.11533788000 1 .6286006 lOE-04 
0.00337362315 0.03993516120 -394.28806800000 



a_d = 

-0.09843658160 -1.445114440E-06 3.514463130E-05 
-3.034783470E-05 -2.06698853000 -0.00511193285 
0.03527477130 -0.24432330100 -92.26523090000 



a_e = 

-0.00451416806 

0.00000000000 

0.00000000000 



0.00000000000 
-0.93336550500 
0.00000000000 



0.00000000000 
0.00000000000 
-15.61232640000 



5.5 Minimal parameter inputs 

The minimum requirements for an input file that contains the MSSM couplings and mass pa- 
rameters are as follows. 

• Gauge and Yukawa couphngs must be specified in the SLHA blocks GAUGE, YU, YD, and 
YE. 

• Quark and lepton mixings must be specified in the the SLHA blocks VCKM and UPMNS. 

• Soft mass parameters must be specified in MSQ2 , MSU2, MSD2, MSL2, MSE2, TU, TD, TE 
and MSOFT. Only non-zero elements of the mass matrices need be specified. For MSOFT, only 
the gaugino soft mass parameters Mi, M2, and M3 must be entered. 

• Parameters jj,, tan 13, and the Higgs vev v are to be entered in the HMIX block. The Higgs 
vcv may be omitted. If so, it is approximated by the value obtained by using the current 
PDG value of the Fermi constant Gp. 

All parameters must be specified at the same scale, using the SLHA Q= convention. An example 
file minimal. spc, provided in the distribution, is shown below. 



# Minimal parameter Input template for Super LFV 
BLOCK gauge Q= 1 . OOOOOOOOE+03 

1 3.62831933E-01 # g_Y 

2 6.45857032E-01 # g_2 
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3 1.07837289E+00 # g_3 

BLOCK Yu q= l.OOOOOOOOE+03 

1 1 8.49362382E-06 # y_u 

2 2 3.59564127E-03 # y_c 

3 3 8.70775686E-01 # y_t 
BLOCK Yd Q= l.OOOOOOOOE+03 

1 1 1. 364646 19E-04 # y_d 

2 2 2 . 86579735E-03 # y_s 

3 3 1.36927282E-01 # y_b 
BLOCK Ye Q= l.OOOOOOOOE+03 

1 1 2 . 98272756E-05 # y_e 

2 2 6.16732139E-03 # y_mu 

3 3 1. 037648 17E-01 # y.tau 
BLOCK VCKM q= 1 . OOOOOOOOE+03 

1 1 9.74192032E-01 
1 2 2.25694592E-01 

1 3 3.44018712E-03 

2 1 -2.25645361E-01 
2 2 9.73359257E-01 

2 3 4.06930840E-02 

3 1 5.83567101E-03 
3 2 -4.04191404E-02 
3 3 9.99165771E-01 

BLOCK UPMWS Q= 1 . OOOOOOOOE+03 
1 1 l.OOOOOOOOE+00 

1 2 7.68699045E-08 

2 1 -7 . 68699045E-08 

2 2 l.OOOOOOOOE+00 

3 3 l.OOOOOOOOE+00 
BLOCK Tu Q= l.OOOOOOOOE+03 

1 1 -4.99689032E-03 
12 2 . 26225278E-09 

1 3 3.24991717E-08 

2 1 9 . 57689089E-07 
2 2 -2. 11533788E+00 

2 3 1.62860061E-04 

3 1 3.37362315E-03 
3 2 3.99351612E-02 
3 3 -3.94288068E+02 

BLOCK Td q= l.OOOOOOOOE+03 
1 1 -9.84365816E-02 
1 2 -1.44511444E-06 

1 3 3.51446313E-05 

2 1 -3 . 03478347E-05 
2 2 -2.06698853E+00 

2 3 -5.11193285E-03 

3 1 3.52747713E-02 
3 2 -2.44323301E-01 
3 3 -9 . 22652309E+01 

BLOCK Te Q= 1 . OOOOOOOOE+03 

1 1 -4.51416806E-03 

2 2 -9.33365505E-01 

3 3 -1.56123264E+01 
BLOCK MSOFT Q= 1 . OOOOOOOOE+03 

1 1 . 03432998E+02 # M_l 

2 1.93115207E+02 # M_2 

3 5.68025363E+02 # M_3 
BLOCK MSL2 Q= 1 . OOOOOOOOE+03 

1 1 6.16982264E+04 

1 2 l.OOOOOOOOE+00 # The source of LFV lies here. 

2 1 l.OOOOOOOOE+00 

2 2 6. 16872834E+04 

3 3 5.86013300E+04 
BLOCK MSE2 Q= 1 . OOOOOOOOE+03 

1 1 4.20218854E+05 

2 2 4.20196836E+05 

3 3 4.13987017E+05 
BLOCK MSQ2 Q= 1 . OOOOOOOOE+03 

1 1 7.20764254E+05 
1 2 4.63631320E+01 
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1 


3 


-1, 


. 12707907E+03 


2 


1 


4, 


.63631320E+01 


2 


2 


7, 


.20449561E+05 


2 


3 


7, 


.80647919E+03 


3 


1 


-1, 


, 12707907E+03 


3 


2 


7, 


.80647919E+03 


3 


3 


5, 


.22563390E+05 


BLOCK 


MSU2 


5= 


l.OOOOOOOOE+03 


1 


1 


7, 


. 68519383E+05 


1 


2 


1, 


.65077717E-09 


1 


3 


5, 


.78795005E-06 


2 


1 


1, 


.65077717E-09 


2 


2 


7, 


. 68512439E+05 


2 


3 


2, 


. 90046080E-02 


3 


1 


5, 


.78795005E-06 


3 


2 


2, 


.90046080E-02 


3 


3 


3, 


.76076532E+05 


BLOCK 


MSD2 


Q= 


1 . OOOOOOOOE+03 


1 


1 


6, 


.90390410E+05 


1 


2 


-7, 


.61272825E-07 


1 


3 


8, 


.84248145E-04 


2 


1 


-7, 


.61272825E-07 


2 


2 


6, 


.90385750E+05 


2 


3 


-1, 


.28617580E-01 


3 


1 


8, 


.84248145E-04 


3 


2 


-1, 


.28617580E-01 


3 


3 


6, 


.79889912E+05 


BLOCK 


HMIX 


Q= 


l.OOOOOOOOE+03 



1 3. 998239 10E+02 # mu 

2 l.OOOOOOOOE+01 # tan(beta)(q) 

3 2.43605877E+02 # vCQ) 



6 Summary and outlook 

Presently we are in an era in which MEG, a dcidicated — > 67 search, is currently acquiring 
data, along with an operational LHC, an unprecedented variety of future non-collider experi- 
ments, and future B-factories. Accordingly, a need for rapid model-discrimination arises. We 
have introduced SuperLFV , a new SLHA tool for calculating low-energy lepton flavor violating 
observables in the context of supersymmetric models. SuperLFV offers a few unique features, 
including (1) independence from existing SLHA spectrum calculators, (2) the ability to accept 
an SLHA spectrum file generated by an existing specturm calculator, personal code, or by hand; 
and (3) an ability to report observables and couplings as a sum of contributions. Features (1) 
and (2) allow arbitrary model selection. Feature (3) promotes analytical insights. 

SuperLFV is intended to be a well-supported addition to the SLHA library. The initial offering 
is minimally comprehensive. Future enhancements will be demand-driven and may include 
greater accuracy {e.g., wavefunction renormalization of the neutralino and chargino fields, more 
sophistication in handling the nuclear aspects of coherent muon conversion, etc.), the inclusion 
of neglected amplitudes {e.g., Higgs- mediated LFV), more observables, or possibly expansion 
beyond the MSSM. 
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The SuperLFV uses public code from the Apache Gommons library for some mathematical 
operations. 

A MSSM couplings 

In the following, sfcrmion indices are denoted as /. fcrmion indices are j, and neutralino 
(chargino) indices are A. In 4-component Dirac notation, the interactions of a sfermion-fermion- 



neutralino are summarized below. 

-J^uuxo = uiu,{n'i^^PL + n'i^^PR)x°A+i^-C- (43) 

-£,-,^0 ^ djd,in''j-^PL + n'}f^Pn)^"A + h.c. (44) 

-C,,^o = g/e,(nJ^-^Pi+n^«^P;j)x!\+h.c. (45) 

-Cp^^o = y,DLj{ntjAP^+<UPR)^A + '^-C- (46) 

^ V2gyY^N:^l}'+^^* + y^N\^t (47) 

^ V2gyYQNl,Ul^* +V2g^tl^N\^Ui=* +yrNAAUi'-'+^^ (48) 

n^^^ ^ X^5FljiVli[/^{(^+^)*+y^%l3[/f (49) 

n% ^ ^^5FFQ^^llt/f* + ^/2524^l2t^^+yf ^^A3t/f+'^* (50) 

^ V2gyY,N\,ul^'+'^* +y:^N-\,Ut (51) 

^ V^5yi^LiViiC/^^* + V2ff<iVl2t//^*+yrAr^3C/f+'^* (52) 

= (53) 

nl-^ ^ V2gYYr.NX,Ur + V2g2tl^NX,Ur (54) 

The sfermion-fermion-chargino interactions are as follows. 

-^dux ^ rf/%(c^jk^L + c^«4Pfl)xi+h.c. (55) 

-Cudic = uidjic'i^^PL + c'jjMxA + '^-c- (56) 

-Uex = i>iej{4}APL + 4j'MXA+^-c- (57) 

-£e.x = eiULMUPL+c-ifAPR)x\+^-c. (58) 

c/j^A = -yjfVi^f/i'* (59) 

c]'j!4 = 92UAiU^'*-yrUA2U^^'+'> (60) 

C/^A ^ -yf C^l2t/i** (61) 

c^«4 = 52V^iC/-^^*-y;^t/l,f/-^(^+^)* (62) 

clU = -ye'^*C^l2C^f* (63) 

<j!4 = 52V:4iC/^'* (64) 

= (65) 

clfA = Q^UaiU'J* -vTUA^ui^'^^'^* (66) 
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From here, one may calculate the effective couplings. The electromagnetic form factors 
of eqn. (21) are generated at one-loop order by neutralinos and charginos. For muon decay, 
these operators are to be evaluated at the muon mass scale, i.e. A^2R{'f' - 
electromagnetic form factors couplings are the following. 



Following [5], 



^nl^'ijA 2-9x + 18x^ -llx^+ 6x^ \nx 



1^ 5767r2 ml - 



^1 n — ^1 1 



1 CjfAC/M 16 - 45y + 36y2 - ly^ + 6(2 - 3y) In y 
5767r2 m? (1 - yY 



A 



IR — ^il\L^R 

1 '^TiA'Ti'TjA 1 - 6x + + 22-3 _ ^ 



u 

2^ 327r2 ml 6(1 - a;)^ 

^ 1 "7g4"/M^x"^ 1 - x2 + 2a: In X 
327r2 m^^ml (1 — a;)^ 

1 c^t>JM 2 + 3y - 6;/2 + ^3 ^ y 

327r2 m?^ 6(1 - y)^ 

1 ^ul^'ijA'^x-^ (-3 + 4y-y^~2liiy) 

327r2 m^-ml (1 — y)^ 



where 



(67) 



(69) 



- AI'lIl^r (70) 

Above, x = xiA = TO|o^/m|^ and y = yjA^ TO|(^/m|^. 
The effective couplings for Z-exchange are 

F'l' = :-]-^n'jf2nnBiNA3NB3-NAiNB4){FiAB + 2GiAB) 

-^2^iTa^T,b (^Iva2Vb2Fiab - Ua2Ub2GiaB^ (71) 

= ^n'j^Xn'j^jBiNA3NB3-NA4NBi){FiAB+2GiAB) (72) 



FjAB = lnxjA + ^^(^^^^^^^^] (73) 

XI A - XiB \ I - XI A 1 - XiB 

^nd GiAB = ' / x,^lnx,^ _ a:zBlnx,^ ^^^^ 



Here, m? is when used in the neutralino contribution (marked by uha couplings) or 
when used in the chargino contribution (marked by cha couplings). 

For the box diagrams of — >■ 3ej, the effective couplings are as follows [5]. 



^1, 



2"/fA"jM''./iB''/jB"^4(Ai3/j) +n/fA?^,/;/A'ijiB%"s ' m^'X'm^%h(ABij) 
1 

2' 



-^'^?fA"jjA<jB"?iB ■ m^\m^lh{ABij} 

+ 4C/fACjjACj;;sC/Js'/4(AS/j) - ^Cjf^cf^^cf^gCjfg ■ m^-m^-J^(^ABij) (76) 
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+C?fACj;;-^4«BC/jS • ^"xb ^4(Ai3/J) (77) 

e^^li = \nulnjfAn''j)j^n]'.j^-m^o^m^oJ^ABiJ) (78) 

B'aR = ^aiU^fl. where a e {1,2, 3, 4} (79) 

The loop integrals /4 and J4 are 

r d-^fc 1 

ihiABij) - J (fc2 _ m?^)(fc2 - m?^)(fc2 _ m?p(fc2 - m? ) ' ^ ^ 

^^^-^""^■'^ 7 (27r)4(fc2-m|^)(fc2-m|^)(fc2-m?)(fc2-m2^)- ^^^^ 

The box diagrams of muon conversion have the following effective couplings [5]. 

D""- ^ Dl^ + Dl"- (82) 

e -'-'n - ^\^I2A^IlB^JqA^JqB ^ I2A^ IIB^ JqA^ JqB l'^ i{ABIJ) 

-|('^/2A%?i3«V49B ~ "7?A"/lS^4gA"j?s) ' m^<Xrn^%h(ABI J) (83) 

= ^^^k^fl where a e {1,2,3,4} (84) 
e^^c"- = -gC^fACnBCjiACjtB'^4(Ai3/j) + ^c^flcJfsCjflc^^fB •Wx^'Ti-^ (85) 
e^^c"- = Ic'j^^14?b4ia4TbMabij) - ^c'j^X4^b4\a4Tb ■ (86) 
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